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The rate of the nuclear reaction 3 He + p~ —> 3 H + 7 + has been calculated using both the 
elementary particle model (EPM) approach and the impulse approximation (IA) approach. Using 
the EPM approach, the exclusive statistical radiative muon capture (RMC) rate for photon energy 
greater than 57 MeV is found to be 0.245 s~ and the ordinary muon capture (OMC) rate to be 
1503 s _1 . The IA calculation exhibits a slight dependence on the type of trinucleon wave functions 
used. The difference between the IA and EPM calculation is larger for RMC than for OMC. To 
resolve the difference between the two approaches a more detailed investigation including meson 
exchange corrections will be required. 
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I. INTRODUCTION 


A recent TRIUMF experiment [lj,|2j] designed to measure the rate of the radiative muon capture (RMC) reaction 
3 He + p~ —> 3 H + 7 Y Vfj, has sparked a renewed interest in this process. Since it is more sensitive to the nucleon 
pseudoscalar form factor gp than its non-radiative counterpart (Refs. |j]Q] are two reviews on ordinary muon capture), 
it is an ideal candidate for checking the value of this form factor which is theoretically predicted by PCAC (Partial 
Conservation of Axial Current). With the experiment on-going, it is necessary to have a modern theoretical calculation 
of the process to interpret the anticipated experimental results. In this paper, we have calculated the rate of the process 
using two perspectives: a nuclear perspective via the elementary particle model (EPM) and a nucleon perspective via 
the impulse approximation (IA). 

The only similar calculation for this process was done by Klieb and Rood about twenty years ago but the 
accuracy of their calculation is constrained by the facts that 1) the trinucleon wave function they used is inadequate by 
today’s standards, 2) some of the nucleon momentum terms were handled in an approximate way and 3) they did not 
use the full Adler and Dothan amplitude but included only some of the Adler and Dothan terms. As a consequence 
of better computer technology and better methods of calculation of realistic trinucleon wave functions, it has been 
possible in this calculation to improve on the approximations they made. The analogous nonradiative, or ordinary 
muon capture (OMC), reaction has been considered by several authors. See for example Refs. |?] [nj, and references 
cited therein. 

In section [n| we will briefly discuss various hypotheses governing the weak hadronic current and the Adler and 
Dothan 0 procedure which provides terms in addition to terms generated by a naive insertion of photons on 
each particle with charge or magnetic moment (see also Ref. [[l2| for early work on this topic and Ref. [Q for an 
alternative presentation of the Adler and Dothan procedure). The elementary particle model approach is discussed 
in section III, where both 3 He and 3 H nuclei are treated as single entities with internal structure revealed only by 
the phenomenological nuclear form factors taken from experiments. We will then discuss the impulse approximation 
approach in section IV. The essence of the impulse approximation is that it regards the radiative capture process 
as taking place on the constituent nucleons. Assuming that the nucleons are free, one then uses a trinucleon wave 
function to integrate out the internal degrees of freedom. The resulting amplitude will be one which only depends on 
the EPM, or external, degrees of freedom. The results will be presented in section |v| and a summary in section VI. 
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II. THE TRANSITION AMPLITUDE 


The fundamental terms of the transition amplitude are obtained by inserting a photon on all particles carrying 
charge or magnetic moment. Figure ( 0 ) shows the Feynman diagrams corresponding to these terms. Note that this 
general approach is common to both EPM and IA so that the form of the transition amplitude derived here will be 
applied to the nucleus in the EPM, but to the nucleon in the IA. 





FIG. 1. The fundamental radiating diagrams 

The amplitudes corresponding to the diagrams in Fig. (Q) are 

Ml = u{v) 7 a(1 - 'y 5 )u(n)u(Pf)W a (Q H )S F (P i - k)Q l u(P l ) 

M -2 = u(v ) 7q (1 - 7 5 )u(fi)u(P f )Q f S F {P f + k)W a (Q H )u{Pi ) 

M 3 = M(»7a(l - l 5 )u{n)u{Pf){ m2 _ _ k ^ 2 {-i){2Q H - k) ■ eGf }^^7 5 M(P t ) 

M 4 = uiy ) 7a (l - 7 5 )Mm - k){-iflu(n)u{P f )W a (Q L )u(Pj (1) 

where S F is the Feynman propagator for spin i particles ( S F (Pi — k) = for example), v, /i and k are the 

four-momenta of the neutrino, muon and the emitted photon respectively and W a (Q) is the weak hadronic vertex 
which is parameterized by four form factors 

W a (Q) = G vl a + G M i<J af + G^V + G P 7 5 — (2) 

2 M n m 

cr a0 = % -( 7 V- 7 V) (3) 

with all the G;’s functions of Q 2 , the square of the momentum transfer at the weak hadronic vertex. Specifically, 
we denote Q H = /i — v as the momentum transfer at the hadronic vertex when one of the hadrons is radiating and 
Q L = fi — v — k as the momentum transfer when the lepton is radiating. Qi(f) = ie^f^P ^p-a Xp k p e\ where e^), «*(/) 
denote the electric charge and anomalous magnetic moment of the initial (final) particle. The induced pseudoscalar 
coupling, Gp, which originates from the pion pole term has the form |lj|| 
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G P (Q 2 ) = 
£ = 


2mM n G a{Q 2 )(X T e ) 


m 2 - Q 2 

ml f G„(Q 2 )/G„(0) 1 
-Q 2t G A (Q 2 )/G A (0)< 


( 4 ) 


which comes from PCAC and the Goldberger-Treiman relation for ordinary muon capture. The quantity e can be 
regarded as a constant over the Q 2 concerned [^] |j. The “PCAC” value of Gp is defined as e = 0. G is the pion-i-f 
coupling constant. In the EPM, the hadronic vertex is at the nuclear level. Therefore, Pi = P3 He (four momentum 
of the 3 He nucleus) and Pf = P3jj (four momentum of the 3 H nucleus). In the IA, the hadronic weak interaction 
operator acts on the constituent nucleons and so Pi = p (four momentum of the proton) and Pf = n (four momentum 
of the neutron). Note also that M n is the mass of nucleus in the EPM but the mass of nucleon in the IA and m (m n ) 
is used to denote the mass of muon (charged pion) throughout. 

However, the sum of these four diagrams is not gauge invariant (GI) and does not satisfy CVC and PCAC by itself. 
Extra terms must be added in order for the whole transition amplitude to be gauge invariant and to satisfy CVC and 
PCAC up to a desired order. The Adler and Dothan procedure is used to generate these terms up to O(k 0 ) via the 
GI requirement and O(Q 0 ) via the CVC and PCAC hypotheses. The extra piece of amplitude that is required is 


AM 

AMi 

A M 2 
am 3 


AMi + A M 2 + A M 3 
-u{v) 7„(1 - 7 5 )»M«(^){Gm’’' t “ C — 


ia — t — + Gp— 7 5 
2 M n m 


}«(^<) 


r 2 k a 'i 

("haO- - 7 5 MvMPf){ - 2(G^7m + G' A l^){k a e^ - Fe“) - G' v (K f - K^-^-ia^k^juiPi) 

M7a(l - 7 5 )u(n)u(P f ){ - 2 (G'y^ + G' A ^l 5 )g^Q H ■ e - 2mM f^ + £ ) 2 Q" ' 
t to 2 - Q l m 

■G' M icJ apC ^f{2Q H ■ e))u{Pi). 


a. 5 
- 7 


( 5 ) 

( 6 ) 

( 7 ) 


( 8 ) 


A prime on a form factor denotes the derivative with respect to Q 2 . Since all the form factors (except Gp which 
contains a pole term) are almost linear in Q 2 over the range of Q 2 concerned in both the EPM and IA, it does not 
matter with respect to which Q 2 the derivative is taken. The pion pole terms of Gp are treated exactly. 

The term AM\ in AM arises when one does a minimal substitution on the hadronic vertex, assuming constant form 
factors, and corresponds to the usual fifth diagram included in previous RMC calculations (Fig.(||)). If the nuclei or 
nucleons were elementary particles with no form factors, this would be the only term needed in AM in order to ensure 
GI, CVC and PCAC of the amplitude. The other terms in AM involving derivatives of form factors are the terms 
to account for the composite nature of both the nuclei (in the EPM) and the nucleons (in the IA). The term AM 2 
is demanded by the CVC and PCAC while the terms in AM 3 are demanded by the gauge invariance requirement in 
the case when form factors are included []. 



FIG. 2. The “minimal coupling” diagram 


1 These are almost all the terms considered by Adler and Dothan 0 except that we have not considered the non-Born terms 
in the pion photoproduction amplitude which are thought to be small. Klieb and Rood |^j6| did not have the last term in AM 2 
and the last two terms in AM 3 . 
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The full amplitude M = Ei=i -4 Mi + AM satisfies GI, CVC and PCAC up to, but not including, terms of O(kQ) 
and this is the amplitude that we will use for later calculations. 


III. THE ELEMENTARY PARTICLE MODEL 

The element ary par ticle model is probably the simplest method to calculate the RMC rate. It was first used by Kim 
and Primakoff |15], HI in calculating the beta decay of complex nuclei and was subsequently used by Fearing |17j and 
Klieb and Rood^H in their RMC calculations. In this approach, both the 3 He and 3 H are treated as “elementary 
particles” of spin p and isospin | up to a small isospin breaking. The details of the structure of the nuclei are 
encapsulated in the nuclear form factors which are determined by experiments done on these nuclei. One of the major 
reasons the EPM calculation is easier than the IA is that there are many fewer degrees of freedom which one has to 
take care of in the EPM than in the IA. The degrees of freedom in the EPM approach are the four-momenta of 3 He 
(Ps He)> 3 H (P 3 h)i photon (fc), neutrino ( v ) and muon (pi) together with their respective spins. The four-momentum 
conservation relation is 


Pan + V + k = Ps He + P- 

The differential capture rate (photon spectrum) is given by 


dT 

dk 


E 

photon polarization 


d{PhaS j k SPaCe) Tr(pM(P 3H: Pane, A'/t) t M(Pa H , Pane, M t )) 


( 9 ) 


( 10 ) 


where d(phase space ) is the differential phase space factor 

ilV ud | 2 G| d 3 Pa H d?v 1 d 3 k 


d(phase space) = C|(0) | ^ 


(27t) 3 (27t) 3 k (27r) 3 


(27t) 4 <5 ( ' 4 ' > (/U, + Pi He ^ Pa H ~ k — v) 


1 |V ud | 2 Gp 2 P 3 ° h 


= C|^(0)| 2 — 1 v ud ' VJp —--S|2fc(l — cos (6)) — 2(m + M t )\ 1 f 2 (k, 9)di>k 2 dkdk 
2n° 2 k 

cos (9) = v ■ k 

,, 2 k(m + M t ) - m 2 - 2 mM t 

v(k 1 9) = 


2k(l — cos (0)) — 2 (to + M t ) 


Mt = 2808.66 MeV. 


( 11 ) 

( 12 ) 

(13) 

(14) 


Here </> At (0) denotes the muon wave function at the origin, C = 0.9788 is the correction factor that accounts for 
the non-pointlike nature of the nucleus, and V ud is the CKM matrix element which connects the up and down quark, 
with V ud = 0.9735 ±0.0008 [[To). Gp is the Fermi coupling constant, p is the density matrix which describes the initial 
spin configuration of the muonic atom, and M(Ps H , Pajje, Mt) is the transition amplitude M that was discussed in 
section (jnj with the following name changes 

«(Pi) -> u ( Pa He) 

u(Pf) -> 

M n —> M t on the hadronic vertex 
Gi —> Fi 


(15) 


(16) 


where Fj,.i = V. M, A are the nuclear form factors which are parameterized as 

Fi = F t ( 0)(1 + pP 2 Q 2 ) 
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with Fy( 0) = 1, Ry = 1-94 fm [R[, Pm( 0) = K 3 He — K3 h = —8.369 n.m — 7.913 n.m, Rm = 1-72 fm Q, Fa(0) = 
1.212 ± 0.004 [g, R a = 1.703 fm0. 

The photon polarization vector e\. which is defined by e\ = -^(x + iXy ), with k = z, has the property 


k x e\ = —iXke \• 

Thus the Ephoton polarization in E( h © can be replaced by Ex=-1,+1- 


(17) 
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The big advantage of the EPM is its simplicity. It also includes, to some extent at least, part of the meson exchange 
corrections which are missing in the IA. It however does have a major flaw as applied to a two step process such as 
RMC in that the intermediate states have to be treated as elementary particles as well. Thus effects coming from 
excitation of the intermediate nucleus, which are implicitly partially included in the IA, are not included in the EPM. 
A full investigation of this is beyond the scope of this paper, though see Ref. (l 8 j . However one should keep this in 
mind as a caveat with respect to the EPM. 


IV. THE IMPULSE APPROXIMATION 


The impulse approximation method provides a simple “microscopic” picture of the nuclear reaction in terms of 
nucleons. In this picture, the constituent nucleons inside the nucleus are approximated as free (this is probably a 
good approximation as the binding energy of the trinucleon system is ~8 MeV, which is about 0.3% of the mass 
of the nucleus) and the nuclear reaction 3 He + fi~ —> 3 H + + 7 is viewed as the sum of its nucleon counterparts 

p + —*■ n + + 7 (that is, only one-body currents are considered and two-body meson exchange currents are 

neglected or put in later as a correction). The extra degrees of freedom which arise from considering each nucleon 
of the nucleus instead of treating the nucleus as a whole are integrated out using realistic tri-nucleon wave functions. 
Given these assumptions, the whole problem boils down to separating the EPM and non-EPM degrees of freedom 
and finding the IA equivalent (denote it as M ia ) of M(Psn,P3n e , M t ) in the EPM. More explicitly, we want the 
M(P 3H , P>He, M t ) in Eq. (10) to be replaced with Mi a in order to find the IA version of The relationship between 
M(P 3 H ,P 3 He , Mt) and M ia is: 


M{P3 H ,Ps He ,M t ) <-> 3 j {{2Trf6 {3 \p' a -p a )(2Tr) 3 S {3) (^ a -q a + ^{u + k-jl)}ipS n (p' a ,q , a )ip3 H e(Pa,q a ) x 

d 3 q^ d 3 f a d 3 q a d 3 Pa 


M(k' a ,k a , M p ) 


(27t) 3 (27t) 3 (27t) 3 (27t) 3 


} = Mu 


(18) 


where k a ,kp,kry (k' a ,k'p,kly) denote the four momenta of the three initial (final) nucleons a, (3 and 7 . The (three) 
momentum transformation separating the EPM and non-EPM degrees of freedom is 

P = ka T ky 3 A kry 


<?« = \k a - l -kp - 

Pa = — —kj 


(19) 


for the initial nucleus, with an identical transformation law for the final nucleus. It is obvious that P (a EPM 
degree of freedom) is the center of mass momentum vector of the nucleus and q a (a non-EPM degree of freedom) 
is the momentum of nucleon a (spectator) with respect to the center of mass momentum of the other two nucleons 
(subsystem) while p a (a non-EPM degree of freedom) is the momentum of particle (3 with respect to particle 7^. 
The 3 comes from the antisymmetrization of the wave function and it allows one to let the current operator act on a 
particular nucleon (chosen to be nucleon a) three times instead of acting on each nucleon of the nucleus. 

The parameterization of the nucleon form factors is exactly the same as that of the nuclear ones (see Eq. ©)■ 
For convenience in notation, we change Fi —> g t and Ri —> r, in Eq. (16) to denote the nucleon case. The various 
parameters for nucleon form factors aregy (0) = 1, ry = 0.7589 fm 0. 9m{ 0) = Kp-i t„, r M = 0.8781 fm |J, £u(0) = 
— 1.267 ± 0.0035 Jl|], r A = 0.6580 fm §. 

The momentum space trinucleon wave functions |2(i| , pl|] are realistic wave functions derived from the Faddeev 
equation (see, for example, Ref. (^2|) with different model potentials. Each one of them has 22 channels which contain 
all possible states up to and including J = 2, where J is the total angular momentum of the subsystem particles. 
They can be written as, 


I = ©V%(Pa,<7a)l ic) | P) 

i c 


( 20 ) 


2 p a and q a will sometimes be denoted as p and q when no confusion arises. 
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where i c is the channel number and | P) is the center of mass momentum of the trinucleon system which we will not 
write explicitly from now on. The coupling scheme of the channels is, 


| 4) = | ((L a l a )C a ,(S a s a )S a )J)\ (. I a i a )l) 
= \i c (J))\i c (I)). 


( 21 ) 


The spin and angular momentum part is given by | (( L a l a )C a , (S a s a )S a )J) with L a (S a ) the angular momentum 
(spin) of subsystem (/I 7 ) and l a (s a ) the angular momentum (spin) of particle a (the spectator particle). These are 
coupled to form C a ( S a ) and then to J = \. The | ( I a i a )T) is the isospin part, with I a being the subsystem isospin 
which couples with i a , the spectator isospin, to form X = 

It is now clear that ip3n e (p a , q a ) in Eq. ( jig ) is just the (p,q\ projection of Eq. (|2c| ) 


e(p a ,q a ) = ^2^i c {p a ,q a ){p,q | ic) 

i c 


where 

(p,q\i c ) = ( 'y£i(p . 0) ® X%}) 


( 22 ) 


(23) 


Note that Eq. ( |26| ) below defines the bipolar harmonic Y£ l (p,q). \ and V i R the above equation are spinors and 
isospinors respectively. Their coupling method is exactly the same as the bipolar harmonic in Eq. (p 6 |). 

To put M ia in Eq- @) into a useful format, we have to expand M(k' a , k a ,M p ) non-relativistically in powers of the 
struck nucleon momentum k a (which equals q a upon setting the initial center of mass momentum of the trinucleon 
zero) and the <5 functions into angular momentum eigenstates. Upon setting jl = 0 and denoting s = v + k the 5 
functions can be expanded as (see Ref. |2^] for definitions of spherical Bessel functions ji(x)), 


(27r) 3 ^ (3) (p' a -p a ) 

(27r) 3 <5 {3) (<f a - q a + |s) 


\rh m 3 
li I 2 


(x,y) 


(27T) 3 


S(.P' a -P<x) 


P 


2 

O' 


^(-i) i V2nrLY?t°tf a ,pa) 


(24) 


'(4tt) 5 (24 + 1)(24 + 1) ( l 3 


^ V ^ 

li,m 3 1 

i—1—3 

X J jiM'cc r )jiMocr)ji 3 {^sr)r 2 dr 
E ( ; | Y™'(x)Y™(y) 


h fo 
0 0 




(25) 

(26) 


The notations used here are the same as Brink and Satchler 124] except that the Clebsch-Gordon coefficients are 


denoted by 


J 

M 


Ji Ji 

M i M 2 


as opposed to (JM \ J\ J 2 MiM 2 ). The next thing after the expansion is to couple all 


the spin and angular momentum operators in Eq. ( |l 8 | ) into tensors of rank 0 or 1. Since the total angular momentum 
of both initial and final states is i, there is no need to couple the operators into tensor of other ranks. There is also no 
need to couple operators into odd parity quantities as both the initial and final states are of even parity. Recognizing 
the total angular momentum of the trinucleon system in the IA as the spin in the EPM, one can easily see that any 
operators of rank 1 in the IA correspond to (within a factor) a matrices in the EPM and operators of rank 0 in the 
IA correspond to identity hadronic operators in the EPM. 

All the coefficients in M(k ' a , k a , M p ) are expanded to O(jj-) except that of gp. The kinematic endpoints of RMC 
are quite close to the poles of gp and thus might make its value large at those places. Therefore, coefficients of gp 
are expanded to 0((jj-) 2 ). 

The correspondence between the IA and the EPM for all forms of operators up to first order in momentum is shown 
below. Operators of higher order in momentum will not be shown owing to the lack of space. Note that the [.. .]’s are 
actually reduced matrix elements between the initial and final states (i.e. results of integration of “internal” degrees 
of freedom) and the numbers inside denote some specific spin and angular momentum combination. They will be 
defined in Eq. ([l3|). For now, it is sufficient to note that the first digit of [...] is related to the nucleon momentum 
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q (for example 0~ no nucleon momentum, 1~ q) and the second digit comes from the spherical harmonics of the <5 
function. These two terms couple together to an angular momentum value represented by the third digit. The fourth 
digit is related to the hadronic spin matrix and the subscript is the rank of the whole reduced matrix element. 1 (or 
sometimes denoted as lo) is defined as the hadronic identity matrix element in both the IA and the EPM. It differs 
from [ 1 ]° (see Ref. |2l|), although they are related. 


IA after coupling and reexpressing in EPM format 

1 «- [((0,0)0 0 0) o ]l 

a ■ v —-^=[(( 0 , 2)2 0 l)i]<r • sv- s + {-^=[(( 0 , 2)2 0 l)i] + [(( 0 , 0)0 0 l)i]}<r • v 

9 ■ <?- {~\j\ [((1,1)2 ® l)i] - |[((1,1)1 ® l)i] - E[((i,1)0 0 l)i]}<? • a 

q ■ V <-> - — [((1,1)0 0 0)o]u ■ s - -^=[((1,1)1 0 0)i]cr • V X S 
Bxq-v <-► {-y^[((l,l )2 0 l)i] + ^[(( 1 , 1)1 0 l)i] + ^[(( 1 , 1)0 0 l)i]}i? x s • <j + 

q-va-u «-► {y^[((l>l)2® l)l] - ^[((l-1)° ® l)i] - ^ [((1, 3 ) 2 0 l)i ]}B ■ uv ■ s + 

!) 2 ® l)i] + |[(( 1 , 1)1 ® l)i] - \\j\ K(l> 3)2 0 l)i]}(? • vu ■ s + 

{_ ys [((1,1)2 ^ ik] - 1)1 ® i)i] - ^ 3)2 0 i)i]} ^ • *+ 


[((1,3)2 ® l)i](j • su ■ sv - s T z^[((l, l)10l)o]u'i?x j 


(27) 

(28) 

(29) 

(30) 

(31) 


(32) 


Note that while B on the left hand side acts on the spin of the spectator nucleon, B on the right acts on the entire 
nucleus and that u and v are mutually commuting vectors that are not concerned with the internal momenta (i.e. 
not p nor q) and commute with B. Using Delorme’s (2^] notation, [ 1 ]° = [((0, 0)0 0 0)o], [c ?] 0,1 = [((0, 0)0 0 l)i], 
[ct ] 2,1 = —[((0, 2)2 ® l)i], [ct] + = [ct ] 0,1 + v^cr ] 2 ’ 1 and [B]~ = [a ] 0,1 — -^[a] 2,1 . The precise relationship between 

[zP ] 1,1 and the reduced matrix elements defined here is unclear but it has a magnitude of the order of [(( 1 , 1)1 0 0 )i] 
or [(( 1 , 1)1 ® 1 )i]. As one will see later, these two matrix elements are very small. 

The definition of [((a, 6 )c[a] 0 d) e ] (a function of s = ||z? + k ||) is 

[((a, 6 )c[a] 0 d) e ] = 




E E (-l) il z Zl - i 2 +b y'{p 2 dpr 2 dr Jf) (^r){^,(p,g') ^{q'r)q' 2 dq'} 

ll 5^2 -,L 1 


{ipic(P> q)jh(qr)q 2+a dq}} 


h h 
0 0 


F(Zi,J 2 ;a, 6 ,c; i' c ,i c ) 


( 2 Li + 1 )( 2 12 + l)( 2 /i + 1 ) 


26 + 1 


(*c(^)ll{(T i ° 1 , il (p / ,p)0 Yfatf, g))®(l 0 ®T d )} e ||i c (J))(z'(J)||(lo®f')||z c (J)) 


(33) 


where Td(e) = lo for d{e) = 0 and B for d(e) = 1 ; r is the isospin operator. Notice a specifies the (mass) dimension 
of the matrix element. When a is not shown explicitly on a reduced matrix element, d = a; that is [((a, 6 )c® d) e \ = 
[((a, 6 )c[a] 0 d) e ]. F(li, Z 2 ; a, 6 , c; i' c , i c ) is defined as 


F(h,l 2 ',a,b, c; i' c ,i c ) = 




Y^iv, xTY^r (v, xK’Z”^ X)dvdx 
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(34) 


J_(_iyi+k+fc/ l 2 b 
4n y 1 a c l o 




( 2?2 + 1)(2 a + 1 ) 


26 + 1 


Sl^r 


v, x being some dummy angular variables. 

(i(,( l 7)||{(y^i Ll (p',p ) 0 Yi? i(q', q)) 0 (lo <8 Td)} e ||i c (j7')) is the spin and angular momentum part of the reduced 
matrix element between the hclion and triton channels. Its calculation is tedious but standard. 

W)\\{(Y^, Ll ip',P) ® Yfrtf, q)) 0 (1 0 0 T d )} e ||z c (J)) = 


47t(-1) £ ' 



-\-‘2l-\-c-\-L' -\-S r -\-^-\-S r -\-d 


e 

c 

d 


C £ 
V l 



{(2 C + 1)(2£' + 1)(25 + 1)(25' + l)(2c + l)(2e + 1)} 
? | s' } S ss'(l\\T d \\l)(L' 0\\Y° uLl (p',p)\\0 L). 


1 

2 


Note also that 

^(-1) l V2Li + 1 (L' 0\\Y° i!Ll (p',p)\\V L) = ±.5 ll> . 

L i 

(?(,(!)||lo 0 T\\i c (X)) is the isospin contribution of the reduced matrix element, which equals 


2 (~ 1 ) 7 ' 


1/2 1/2 1/2 

1/2 1/2 I' 


5w{\wn\) 


(35) 


(36) 


(37) 


By following through all the procedures mentioned in this section, it is possible to match the IA amplitude piece by 
piece with its EPM counterparts and thus make a direct comparison between each piece. In other words, we have 
arranged the non-zero terms of Mi a to depend only on the EPM degrees of freedom. We then obtain ^ via Eq. (0) 
with M(P 3 H ,P 3 He ,Mt) replaced by M ia . 


V. RESULTS AND DISCUSSION 


A. Results 

We have calculated the rate for OMC and the photon spectrum and integrated rate for RMC in both IA and 
EPM approaches, using the formalism described in the preceding sections. Table (jl!) shows various RMC and OMC 
statistical rates, with the final results given in the last three columns. Figure (^) shows the IA RMC spectra for 
various wave functions together with the EPM spectrum. 

One sees immediately that the IA results are significantly lower than the EPM results both for OMC and RMC. 
This is consistent with the OMC results of Refs. For OMC the difference is understood to arise from various 

meson exchange corrections which are included implicitly in the EPM approach but which must be put into the IA 
by hand as specific corrections. [|l|. Presumably a similar explanation holds for RMC. 


Potentials 

OMC rate (statistical) (s x ) 

rss? (s- 1 ) 


0(*£r); I 

0 (i&); fuU 

0(i &)’! 

0 (£,); full 

k > 5MeV 

k > 57MeV 

Bonn-A 

1368.6 

1368.1 

1367.8 

1357.9 

0.6255 

0.1691 

Bonn-B 

1341.3 

1340.8 

1340.4 

1330.7 

0.6164 

0.1666 

CD-Bonn 

1336.1 

1335.7 

1335.3 

1326.1 

0.6153 

0.1663 

Nijmegen I 

1298.0 

1297.5 

1297.1 

1288.2 

0.6023 

0.1626 

Paris 

1270.9 

1270.4 

1270.1 

1260.3 

0.5932 

0.1602 

AV14 

1271.0 

1270.6 

1270.2 

1260.0 

0.5929 

0.1601 

EPM 

1503 

0.8263 

0.2451 


8 






























TABLE I. Various OMC and RMC statistical rates calculated with gp or Fp at the PCAC value. The numbers in the 
second column are the values obtained by using the prescription up to O(jg-) terms. The numbers in the third column are 

obtained via the correct approach up to 0(^~) terms. The fourth column has values of the OMC rates using prescription 
up to 0((jg-) 2 ) terms. Numbers in the fifth column are the values obtained by the correct approach up to 0((jg-) 2 ) terms; 
the sixth and rightmost columns contain T^t(fc > 5 MeV) and T r s t^t{k > 57 MeV) respectively. 
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FIG. 3. RMC photon spectra from two EPM calculations, one with the full Adler and Dothan amplitude AM and the 
other with only the terms AMi + AA /3 necessary for gauge invariance. Also included are RMC photon spectra from the full IA 
calculation using various model potentials, plus one example with only AAA + AM 3 . All wave functions used have 22 Faddeev 
components and the permutation is projected on the same set of states. The infrared divergent part is not shown. The KR 
values come from Klieb and Rood j^,^| and are their IA calculation results. 
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B. Importance of various ingredients 


We can see from these results the importance of some of the specific ingredients and improvements which we have 
included in the calculations. In particular, table (j|) shows the results obtained for the IA statistical OMC rate using 
four different methods of treating the nucleon momentum operator q. The most common way of treating this operator 
is to replace q with ^ |^(| (second column) and virtually all existing OMC calculations used this method. However, 
this method is strictly correct only for S-waves to first order in nucleon momentum. Using the correct approach 
discussed in section IV and keeping terms up to O(jfc-) decreases the OMC statistical rate by about 0.5 s -1 (third 


column) as compared to the results from the prescription. The smallness of the effect is primarily due to the minute 
contribution of the P-state wave function to the trinucleon wave function There is a difference of about 10 s^ -1 
between using the | approach (fourth column) and the correct approach (fifth column) when terms up to 

are kept. Although in percentage terms it is just about 0.6%, it is much larger than the 0.5 s -1 difference between 
the numbers in the second and third columns. This is perhaps to be expected as the approach cannot be applied 
to terms higher than first order in nucleon momentum. 

Table (|j) and Fig. (|3j) also show the dependence of both the IA OMC and RMC calculations on various trinucleon 
potentials. In general, the Bonn type potentials seem to give higher (and perhaps better) RMC and OMC results than 
the other potentials. To analyze this properly, let us take a look at the three dominant reduced matrix elements, [1] , 
[ct] 0,1 and [tx] 2 ’ 1 , which are shown as a function of s in Figs. (@), (|), and (§). All the curves produced by non-Bonn 
potentials seem to be a bit separated from the curves of the Bonn type potentials, especially in the region when s 
is large. For [1]° the problem may be partly associated with the numerical normalization since [l]°(s = 0) should 
be unity in principle. However, even though one takes this into account (say, for example, by scaling the curves so 
that they agree with each other at s = 0) the values of [1]° for the non-Bonn potentials are still smaller than that of 
the Bonn potentials, as can be easily seen from the fact that the fractional deviation of the reduced matrix element 
among various wave functions is larger at large value of s. 


3 Note that this effect is a genuine effect and not an effect caused by numerical calculation. To prove this, a two channel 
Yamaguchi wave function consisting solely of S-waves is used to gauge the numerical uncertainty in wave function integration. 
There is an increase of the rate (due only to numerical integration) by 0.1 s -1 , out of a total 1600 s _1 , for the correct approach. 
This difference is much smaller than the difference (stemming from errors in both the numerical integration and the - approach) 
of the calculations of the other regular 22-channel wave functions. 
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FIG. 4. Plot of [11° vs. s for different nuclear potentials 
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FIG. 6. Plot of [a] 2 ’ 1 vs. s. 

This seems to have to do with the different binding energies produced by the Bonn and non-Bonn potentials. The 
non-Bonn potentials generally underbind the trinucleon by about 0.5 to 1 MeV. Congleton and Truhlik || pointed 
out that^| [1]° ~ 1 — const(r 2 ) s 2 and that (r 2 ) scales like the inverse of binding energy, and thus argued that one can 
expect a lower value for [1]° when a wave function with a lower binding energy is used For [cr] 0,1 , they further 
argued that since it is a reduced matrix element for one-body currents (this IA calculation contains only one-body 
currents), the Bonn potential’s weak tensor force makes this matrix element large in magnitude. If their analysis is 
correct, this would potentially explain why the Bonn type potentials consistently give higher results in both OMC 
and RMC IA calculations. The quadratic like curve of [o’] 2 ’ 1 is obvious if one notices j 2 {x) ~ fg for x « 1. Since 
[cr] 2,1 is the result of S-D overlap, a tri-nucleon wave function with a larger component of D wave would probably have 
a larger magnitude of [cr] 2,1 which seems to be the case for the generally higher D wave component of the non-Bonn 
potentials (see appendix ([b])). 

Recently Lahiff and Afnan [^7j suggested that the Bonn-type potentials might be a better choice than the Nijmegen 
potential because the energy dependence of propagators is treated exactly (during the evaluation of the potential via 
time-ordered perturbation theory) in the Bonn type potential but the Nijmegen group removes this energy dependence. 

Thus to summarize, the reason that the Bonn type potentials give a higher result for both RMC and OMC may 
primarily be due to its higher binding energy and possibly its weak tensor force. The higher D wave components of 
those non-Bonn potentials may give a boost to [cr] 2,1 but the smallness of [cr] 2 ’ 1 as compared to the other two would 
make its effect on the IA capture rates small. 


4 To see this, expand jo(fsr) in polynomials and note that const > 0 

*To define the term “lower binding energy”, note that appendix (|b|) shows the Bonn potentials have the highest binding 
energy predictions among all the wave functions. 
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The effects of the Adler and Dothan terms can also be seen in Fig. (^]). The upper two curves show the RMC 
spectrum including the full AM of of Eq. (|J) and including only AM\ + AM 3 , the part required by GI. There is 
clearly a significant difference between these two photon spectra in the EPM. The photon spectrum with the full 
Adler and Dothan amplitude is 5-25% higher for photon momentum starting from 20 MeV than the one with terms 
only to ensure GI. 

The figure also shows six IA photon spectra for different trinucleon wave functions using the full set of Adler and 
Dothan terms plus an IA photon spectrum with just terms of AMi + AM 3 necessary to ensure GI. For this latter case 
a trinucleon wave function from the Bonn-A potential was used for the calculation. For the IA case there is almost 
no observable change in the photon spectrum produced by including the full set of Adler-Dothan terms as opposed 
to just those required for GI. 

If one considers IA and EPM spectra with only the terms from A Mi + AM 3 , i.e. only those necessary to ensure GI, 
the ratio of IA to EPM RMC total capture rate is similar to that of the corresponding OMC quantity (~ 83 — 90%). 
However, as a result of increased importance of the full set of Adler and Dothan terms in the EPM case, once these 
additional terms are included, the IA to EPM RMC capture rate ratio drops by more than 10% to around 73%. 

An obvious reason for this increased sensitivity in the EPM case might be the more rapidly varying form factors in 
the nuclear case as opposed to the nucleon case. However there is also a relative sign flip between between Fa{ 0) and 
gA( 0) which might play a role. To investigate this further we have done three things. We have artificially increased 
rv, tm and ta to 1.95 fm (a value comparable to various nuclear radii) and calculated the resulting IA photon spectra 
using these values (Fig. (f7[)). We have calculated another set of photon spectra by flipping the sign of gA in the IA 
(Fig. (||)). We also provide the IA photon spectra when both of these effects are present (Fig. (||)). 



FIG. 7. Effect of increased nucleon radii on the full IA photon spectrum and on the spectrum with only AM\ + AM 3 
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FIG. 9. Two sets of IA photon spectrum using both the increased nucleon radii and the sign change for qa- 


By looking at Figs. ( 0 - 1 ), one concludes that the IA RMC calculation would have shown a sensitivity to the 
additional terms in the full AM if the nucleon radii were of comparable size to the nuclear ones, no matter whether 
the sign of gA is flipped (as in Fig. (||)) or not (as in Fig. (j7|)). The flipping of sign of gA decreases both spectra (full 
IA and calculation with only terms to ensure GI) but the difference between them is relatively small, provided the 
various nucleon radii are not artificially increased (Fig. (||)). 

One notes that there are some additional deficiencies of the Adler and Dothan procedure which may be relevant. 
In particular, O(kQ) terms in AM, which are formally of the same size as some terms which are kept, cannot be 
determined uniquely by GI or CVC and PCAC alone. As a result, some terms in AM involving derivatives of form 
factors are missing. This might be problematic when those derivatives are large, as in the EPM. 

A final remark should be made regarding the comparison of IA and EPM approaches. A major difference probably 
originates in the fact that the IA as described here uses only one body currents. To elucidate this difference one 
needs to perform a more detailed investigation of the interactions of the nucleons using a model, say involving meson 
exchanges and adding two (or more)-body currents. These steps would entail a much more complicated calculation 
than the one here solely considering one-body currents. Congleton and Truhlik |J calculated the MEC contribution 
to the simpler problem of OMC by 3 He and found that IA+MEC prediction of the exclusive OMC statistical rate 
agrees with both the EPM calculation and experiment 28 1.|] Clearly something similar needs to be done for RMC. 
Finally one should note that there are also difficulties with the EPM as applied to RMC relating to the treatment of 
intermediate nuclear states, as was discussed in Ref. [BUI . 


6 They used a trinucleon wave function P9(] that 


was derived in a slightly different way than the ones we use. 
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C. Dependence on Fp 


One of the main motivations for examining RMC is to obtain information about the induced pseudoscalar coupling 
constant and hence we have obtained results for various values of this coupling. Figure 0) shows the EPM calculation 
of the RMC photon spectrum for several different choices of Fp. There is an increase in the total capture rate and a 
slight shift of the peak of the spectrum to higher photon energy as one increases Fp from 0.25 to 1.75 times its PCAC 
value. 



FIG. 10. The RMC photon spectrum dr J^ at calculated in the EPM approach using the full Adler Dothan amplitude with 
various values of Fp, in units of Fp CAC as determined from Eq. (p|) with e = 0. The KR values are those of Klieb and Rood, 
taken from the relativistic calculation of Ref. || which are not shown in Ref. [^| 

Figure ( 0 ) shows the same quantity as Fig. ( 0 ) but for IA calculation using wave functions derived from the Bonn- 
A potential. The qualitative features are essentially the same as for the EPM case, though the absolute magnitude is 
different, as was discussed above. 
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dr/dk [MeV 



FIG. 11. Photon spectra from IA calculations for various values of gp, in units of gp CAC as determined from Eq. (:J) with 
e = 0. Wave functions derived from Bonn-A potential are used for the calculation. 

Figure ( 0 ) shows the sensitivity of the integrated spectrum, i.e. the RMC capture rate (k >57 MeV) with respect 
to variation of Fp (for the EPM) and gp (for the IA). The increase in total capture rate as Fp or gp increases from 
0.25 to 1.75 times its PCAC value is slightly more rapid in the EPM than in the IA. 


19 












FIG. 12. Sensitivity of > 57 MeV) to Fp or gp. The Bonn-A potential was used for the IA calculations and the full 

AM was included for both IA and EPM results. 

Another quantity which is sensitive to Fp in RMC and can in principle be considered, though the experiment is 
difficult, is the photon polarization. This is defined as the rate (or spectrum) for particular photon polarization minus 
that for the reversed polarization divided by the sum. Figure ( |13|) shows this photon polarization P 7 (fc) using the 
EPM for various values of Fp. Clearly there is a very strong dependence on Fp, particularly for the highest energy 
photons. These photon polarizations for different values of Fp all seem to converge to a limit as k —+ 0. This is due 
to the fact that in this limit the amplitude is determined by soft photon theorems. In the usual transverse gauge and 
for initial muon and nucleus at rest, the leading term in the squared amplitude is O(p-) and has the form | | 2 

(from the diagram with the final hadron emitting). It is thus independent of the sign of the photon polarization, 
which makes P 7 (fc) —> 0 as k —> 0. 

Note that, although we have not calculated it explicitly, similar sensitivities to Fp would be expected in the photon 
asymmetry relative to the muon spin. |3C|] 
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FIG. 13. Photon polarization for various values of Fp. 


D. Comparison with other works 

The EPM photon spectrum (Fig. (|To|)) and total RMC rate, given in table (I), are in good agreement with the results 
of Klieb and Rood who obtained a total rate of 0.814 s -1 j(J via a non-relativistically approximated amplitude. 
Thus the extra terms included here but not included by Klieb and Rood seem not to contribute significantly to the 
photon spectrum. The agreement of the EPM OMC rate of 1503 s -1 with Congleton and Fearing’s [|| is also good, 
though this agreement is perhaps not surprising given that the values of the form factors used here, and the basic 
approach, are almost the same as theirs. 

The IA calculation of OMC rate using wave functions derived from Paris or AV14 potentials agrees with that 
of Klieb and Rood. However, the two IA RMC spectra are about 4% lowci^ than theirs. At first it looks a bit 
contradictory that more or less the same IA OMC rates, but slightly lower RMC results, are obtained in this work, 
but a closer look reveals that Klieb and Rood used a lot of approximations evaluating the reduced matrix elements 
for the RMC spectrum which they did not use for OMC. 


7 This already takes into account the fact that they took the value of C = 0.965 while C = 0.9788 is used in this work. 


21 















In particular instead of evaluating the reduced matrix elements directly in terms of s, they expressed s in terms of 
an infinite sum of spherical harmonics of u and k and imposed an artificial cutoff on this expansion. They also did not 
fully square the resulting matrix element. Only products of any two of the most dominant terms and products of one 
dominant and one small term were considered. In expanding the plane wave exp(*s • r), they only included the term 
having jo(z'r)jo(fcr) and they used this approximation as the premise in deriving several relationships between various 

reduced matrix elements for RMC. They also do not seem to have expanded the spinor normalization factors 

(for the neutron spinor) and (for the proton spinor) quite consistently. None of these approximations were 

used here. Also some terms in AM that are present here but not included by Klieb and Rood tend to decrease the 
resulting photon spectrum a bit. Given these differences, the agreement within 4% for the rates is quite satisfactory. 

The EPM photon polarization P 1 (k ) obtained in this work does differ from that of Klieb and Rood quite significantly 
(see Fig. |l3|). The reason behind this is that while the non-relativistically reduced amplitude used by Klieb and Rood 
produces the correct spectrum within a few percent, it cannot produce P 7 (fc) accurately. Fearing |3(J noted that 
the first order contribution to P 7 (fc) actually comes from C>(—S-^) terms in the squared Hamiltonian. Klieb and 

Rood apparently compromised P 7 (fc)’s accuracy by truncating many 0( j^~t) terms when they squared their already 
non-relativistically reduced amplitude. 



VI. SUMMARY AND CONCLUSION 

We have performed a theoretical calculation of the process 3 He + n~ —> 3 H + 7 + using two separate approaches, 
the elementary particle model and the impulse approximation. Our calculation contains a number of improvements 
over the previous ones, namely: 1) The full Adler and Dothan amplitude is used for both the EPM and IA calculations. 
2) Better momentum space wave functions from various nuclear potentials are employed for the IA calculation. 3) 
The non-relativistic reduction of the IA amplitude contains second order nucleon momentum terms for coefficients of 
gp. 4) The nucleon momentum terms in the IA are treated exactly without using the common | approach. 

In general our results agree well, when comparisons can be made, with the older calculations of Klieb and Rood |^||. 
In particular using the EPM model approach, the RMC statistical rate obtained in this work agrees with that of Klieb 
and Rood.The photon polarization P 7 (fc) disagrees significantly with their results, which probably is a consequence of 
the fact that they truncated a lot of terms in the squared amplitude which contribute significantly to P 7 (fc). 

The IA OMC results derived from wave functions of the non-Bonn potentials roughly agree with Klieb and Rood’s 
but the RMC photon spectra from the same wave functions are slightly lower than Klieb and Rood’s. This seems to 
have to do with the fact that they made some IA RMC specific approximations in evaluating their photon spectrum. 

We summarize our results as follows. As expected there is a strong dependence of the results on the value of the 
induced pseudoscalar coupling constant Fp or gp. There is a slight dependence of the IA calculations on nuclear 
potentials. The dependence can possibly be accounted for by the difference in the three body binding energy resulting 
from the different potentials, by details of the nuclear potentials such as stronger/weaker tensor force, etc., and by 
differing partial wave characteristics of the resulting trinucleon wave functions. 

There is quite a significant difference between the EPM and IA RMC calculations. A first look at the spectra 
would suggest that the difference is caused by the fact that the EPM calculation is much more sensitive than the 
IA to the 0(k,k 2 ) terms in AM which are larger in the EPM because of the much larger magnitude of the various 
nuclear radii than their nucleon counterparts. This hints at a poor convergence of these Adler and Dothan terms in 
the EPM and suggests that the higher order pieces which cannot be calculated via the Adler-Dothan procedure might 
be important. Further problems with the EPM have been discussed before in Ref. jL 8 |. Probably the most important 
effect contributing to this difference however are the meson exchange corrections. These accounted for the difference 
between IA and EPM in the OMC case |j|. Such a calculation will be much more complicated for RMC, but is clearly 
needed to fully understand the differences between the IA and EPM approaches. 
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APPENDIX A: M(P 3h ,P s He ,M t ) IN THE EPM 


Using the Dirac representation of the 7 matrices, M(P 3 H , P 3He , M t ) (up to a constant phase factor) can be written 
in the form below which operates on the product space of the leptonic and hadronic spinors. 
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where 07, is the leptonic spin matrix and a is the hadronic spin matrix. N' is the spinor normalization factor which 
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APPENDIX B: WAVE FUNCTION CHARACTERISTICS 


Table (ij) lists several important quantities that depend on the potential used to generate the wave function. They 
are respectively the binding energy given by the wave functions, various partial wave probabilities and the numerical 
normalization (ip \ ip) num of the wave functions. The numerical normalizations of all the wave functions are not unity 
because the antisymmetrization of the wave functions is projected on a finite set of states. The experimental binding 
energy Eb of 3 He is 7.72 MeV and 3 H is 8.48 MeV [[IJ . 


Potential 

E b 

P(S) 

P(S’) 

P(P) 

P(D) 

1 VOn,,™ 

Bonn A 

8.29 

92.59% 

1.23% 

0.030% 

6.14% 

0.994 

Bonn B 

8.10 

91.61% 

1.19% 

0.044% 

7.16% 

0.993 

CD Bonn ( 3 He) 

7.91 

91.61% 

1.35% 

0.041% 

7.01% 

0.993 

CD Bonn ( a H) 

7.96 

91.63% 

1.31% 

0.041% 

7.01% 

0.993 

Nijmegen I 

7.66 

90.31% 

1.29% 

0.065% 

8.34% 

0.990 

Paris 

7.38 

90.11% 

1.40% 

0.069% 

8.42% 

0.988 

AV14 

7.58 

89.86% 

1.15% 

0.082% 

8.90% 

0.987 


TABLE II. Some important quantities of trinucleon wave functions, 
probability of S-wave component of the wave function and so on. 


The binding energy Eb is in MeV. P(S) denotes the 
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